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Abstract 



D ■ 

^ ■ We reanalyse the topological susceptibility assuming the possibility of a significant para- 

magnetic suppression of the three-flavour quark condensate and a correlated enhancement 
CN) i of vacuum fluctuations of ss pairs. Using the framework of resummed xPT, we point out 

^ ■ that simulations performed near the physical point, with a significant mass hierarchy be- 

f— ^ . tween u, d and s dynamical quarks, are not able to disentangle the contributions from the 

lO ! quark condensate and sea ss-pair fluctuations, and that simulations with three light quark 

^— ^ ' masses of the same order are better suited for this purpose. We perform a combined fit 

f^ ■ of recent RBC/UKQCD data on pseudoscalar masses and decay constants as well as the 

i^ \ topological susceptibility, and we reconsider the determination of lattice spacings in our 

framework, working out the consequences on the parameters of the chiral Lagrangian. We 
obtain (S(3; 2GeV))^'^ = 243 it 12 MeV for the three-flavour quark condensate in the chiral 
limit. We notice a significant suppression compared to the two- flavour quark condensate 
S(2; 2GeV)/S(3; 2GeV) = 1.51 it 0.11 and we confirm previous findings of a competition 
^ _• between leading order and next-to-leading order contributions in three-flavour chiral series. 



1 Introduction 

Over the recent years, a joint experimental and theoretical effort has provided new insight concerning 
low-energy dynamics of QCD, through the investigation of vrvr (re)scattering in K^/^ [IHZ] and K ^ 3tt 
decays [8HT2]. It turns out that, in the Nf = 2 chiral limit (m„ = rtiii = 0) [T3j, the breakdown 
of chiral symmetry (responsible for the presence of pseudoscalar Goldstone bosons corresponding to 
pions) is triggered by the quark condensate S(2) = —liiaimu,ma^o{uu). The situation is not so clear 
for the pattern of chiral symmetry breaking in the Nf = 3 chiral limit where m^ also vanishes and 
kaons and r] become also Goldstone bosons [lH[15]. In particular, one would expect chiral order 
parameters (such as the quark condensate) to decrease as the number of massless flavours in the 
theory increases. This suppression is expected to be more and more pronounced as more and more 
light quark masses are sent to zero due to a paramagnetic suppression from the vacuum fluctuations of 
sea-quark pairs [16] . This effect is hinted at by the significant violation of the Zweig rule in the scalar 
sector, and should affect three-flavour dynamics involving kaons and eta mostly, but not (or only 
mildly) the two-flavour dynamics of pions alone [17pi8j. Such a suppression would lead to a decrease 
of the condensate in the Nf = 3 case compared to its value in the Nf = 2 limit. An unambiguous 
test of the pattern of Nf = 3 chiral symmetry breaking comes from lattice simulations, which have 
made significant progress in the recent years, due to many improvements both at the conceptual 
and technical levels. A first avenue consists in computing the spectrum and/or the dynamics of 
light pseudoscalar mesons in order to extract the low-energy constants encoding the pattern of chiral 
symmetry breaking in the two chiral limits (for a recent overview, see ref. |15j). 

Two other alternative strategies are often used in lattice simulations to determine the three- 
flavour quark condensate, both being related to the topological properties of QCD on the lattice. 
The first one consists in studying the infrared end of the spectrum of the Dirac operator directly, 
as the accumulation of its eigenvalues around zero provides information on the value of the quark 
condensate [19p20] . We will not address this type of determination in the present article, but we notice 
that a joint JLQCD and TWQCD lattice study of the distribution of low-lying Dirac eigenvalues 
obtained [21]: S(2)/S(3) = 1.30 ±0.54, where 5](A^j) denotes the quark condensate in the chiral limit 
of Nf massless flavours, in agreement with a paramagnetic suppression of chiral order parameters as 
the number of massless flavours increases. 

A second possibility consists in studying the so-called topological susceptibility, which will be the 
main focus of the present article. One can introduce the winding- number density: 

where a trace over colour indices is performed, and whose integral Q = J dx co can be used to classify 
the gluonic configurations according to topology. In the case of lattice simulations, the winding num- 
ber of a gauge configuration can be determined through the choice of appropriate periodic boundary 
conditions ^19j. Its conjugate variable in the QCD Lagrangian is the so-called vacuum angle 9, 
whose value is strongly constrained by the absence of CP violation observed in strong interactions. 
The value of is experimentally very small {9 < 10^^*^ [22j), which requires either fine tuning or a 
dynamical explanation (e.g., axions |23y24j). 

The topological susceptibility defined as the mean square winding number per unit volume x = 
{Q'^)/V occurs when connecting an observable computed for a fixed topology (i.e., a fixed winding 
number) to the same observable in the 9 vacuum at finite volume [25]. Up to 0{1/V^), this relation 
involves both x Bind the coefficient C4 = [{Q"^) — 3((Q^))^)]/V^ from the 4-point correlator of winding 
number densities (we will come back to this second quantity in a forthcoming article [26j). The 
topological susceptibility is related to the value at zero of the correlation function: 

X(p2) = _i fd'^xe'P-''{0\Toj{x)oj{0)\0) . (2) 



This correlation function can be obtained from the generating functional of QCD by performing two 
derivatives with respect to the local source corresponding to the vacuum angle 9{x). As discussed in 
ref. [27], this correlation function is too singular in QCD for the integral to exist, so that eq. ([2]) is 
an ambiguous notion and has to be renormalised. This problem however does not affect the value of 
X at zero momentum transfer and per unit volume x = x(0)/^) i-6-, the topological susceptibility. 
Indeed, x can be related through Ward identities to correlators requiring no renormalisation: 

x{0) = -'-jdx{0\Taoix)ao{0)\0)-^{0\qmq\0), (J^ = sj'^qil^mq . (3) 

The two terms on the right-hand side of this Ward identity (written in the case of degenerate 
masses among quarks collected in a flavour multiplet q [27]) shows the connection of the topological 
susceptibility with the determination of the quark condensate and the propagation of flavour-singlet 
Goldstone bosons respectively. It is thus of no surprise that this quantity, related to topological 
properties of QCD, is also exploited to determine the quark condensate on the lattice [25ti52] . 

This is reflected at the level of the low-energy effective theory of QCD, i.e. Chiral Perturbation 
Theory (xPT) |13yi4j. One builds a generating functional including the source corresponding to the 
vacuum angle, 6{x). The latter is only affected by axial flavour-singlet transformations, and thus 
connected to that of the (pseudo) scalar source term through: 

[s{x) + ip{x) = so + ipo,0{x) = 9o] ^ [s{x) + ip{x) = VR{so + ipo)Vl,e{x) = eo + argdet(y^yL)] . (4) 

The source 9{x) corresponds to the phase of the matrix collecting Goldstone bosons in the building of 
the low-energy effective theory, yielding Ward identities connecting correlators with (pseudo) scalar 
and winding-number densities. Let us note that a general mass matrix so+ipo can be brought down to 
a diagonal one with real positive eigenvalues rriu, rrid, mg through the axial transformation eq. @, so 
that the effective theory depends on so,po, Oq only through m„, tti^^, rus and 9 = 6*0 + argdet(so + «Po)- 
If one of the quark masses is equal to zero, this transformation can also be used to shift the vacuum 
angle 9 by an arbitrary amount without modifying the theory: the xPT effective potential becomes 
independent of 9, and its second derivative, i.e., X; should vanish in this limit (this can be also seen 
from the Ward identity eq. ([3]) written for several distinct quark masses). 

Eq. dH) yields a natural connection between ^-related observables, the Ua{^) anomaly and the 
rj' meson, connection which can be investigated further by performing a simultaneous expansion in 
momenta, quark masses and l/N^ and promoting the i]' meson to a ninth Goldstone boson [191127] . 
This has led to a detailed study of the topological susceptibility in the pure Yang-Mills theory (see, 
e.g., refs. [33^34]). But even without relying on the large- A'^c expansion, one can determine the low- 
energy behaviour of x(p^), and in particular its value at zero, within three-flavour Chiral Perturbation 
Theory. The leading-order result follows the structure of the vanishing-momentum Ward identity 
eq. ([3]) closely 



X{p') = - E ^S^^ + ^^(3)(m„ + m, + m,) + ip2^0 + O(/), (5) 

p=n0.n Mp -p^ 

where S(3) is (the opposite of) the quark condensate in the Nf = 3 chiral limit |j, (0|a;|P) is the 

o 2 
coupling of flavour-singlet Goldstone bosons to the winding number density and Mp the leading- 
order (LO) contribution to their masses, while Hq is a high-energy counterterm affecting x'(0) but 



^In the following, we will denote S(A^/) and F{Nf) the quark condensate and the pseudoscalar decay constant in 
the chiral limit where the masses of the lighter Nf quarks vanish. According to refs. 13,14 : 

Fo = F(3), Eo = E(3) = F^Bq, F = F(2), E = E(2) = F^B . (6) 



not x(0)- Replacing the first term on the left hand-side by its appropriate expression and setting the 
momentum transfer to zero, one finally gets the widely used and very simple leading-order expression: 

X = S(3)m , 1/m = l/niu + l/m^ + l/m^ , (7) 

with an interesting non- analytic structure going to zero for vanishing iriujTnd or rus- Obviously, an 
accurate determination of the quark condensate from the topological susceptibility requires to go 
beyond a LO analysis. The next-to-leading (NLO) expression for the topological susceptibility was 
recently computed in ref. [35J, by determining the dependence of the xPT potential with respect to 
the vacuum angle (i.e., taking the source term 6{x) as a constant) and performing the corresponding 
partial derivatives. The LO formula ([7]) was exploited by the TWQCD collaboration to extract 
the value of the three-flavour quark condensate |28j from RBC/UKQCD configurations with 2-1-1 
domain-wall fermions [36]: S(3, 2 GeV) = [259(6)(9) MeV]'^ (in the MS scheme). Another study was 
performed by the RBC and UKQCD collaborations based on data at larger volumes using the NLO 
formula [3U], but it was not fully exploited as only a consistency check with the dependence on niu^d 
was performed and the three-flavour quark condensate was not determined by this method. 

In this context, it is interesting to consider the topological susceptibility as a way to extract the 
quark condensate, taking into account the possibility that it is not the order parameter triggering 
three-flavour chiral symmetry breaking. Indeed various analyses of lattice results p ^[371 - li2] suggest 
an overall good agreement between lattice simulations and xPT concerning chiral series obtained as 
an expansion in m,„ and m^i only {Nf = 2 xPT), but significant disagreements concerning the chiral 
expansions in powers of rriu, rud, rng {Nf = 3 xPT). In some cases, small values of the Nf = 3 quark 
condensate and pseudoscalar decay constant fail to saturate the LO chiral expansions. In other cases, 
the convergence of the three-flavour chiral series exhibits sometimes unexpected features, as small 
next-to-leading (NLO) contributions and large next-to-next-to-leading (NNLO) ones, or numerical 
cancellations between LO and NLO terms. 

A scenario where such features could occur was discussed in a series of articles p^HTSlHHtilT] : 
significant vacuum fluctuations of ss pairs, hinted at by the violation of the Zweig rule in the scalar 
sector, could lead to a paramagnetic suppression of LO chiral order parameters (quark condensate 
and pseudoscalar decay constant). In the Nf = 3 chiral expansions of observables such as decay 
constants, masses, form factors. . . these fluctuations would decrease the size of LO contributions and 
enhance NLO contributions numerically (see also refs. |48y49] for a discussion of the r] dynamics in 
the same framework). This would require a particular treatment of chiral expansions, going under 
the name of Resummed Chiral Perturbation Theory (RexPT). It resums higher-order contributions 
in chiral series from the -L4 and Lq NLO low-energy constants encoding the effect of ss pairs on the 
structure of the chiral vacuum and may induce a significant m-s-dependence in the pattern of chiral 
symmetry breaking. It is compatible with the usual treatment of chiral series in the limit where the 
latter are saturated by their LO term, but it allows for a consistent treatment of the series even if 
there is a significant competition of LO and NLO contributions for some of the observables. 

In ref. [17], we illustrated our procedure for energy-dependent quantities, namely form factors, 
such as the electromagnetic form factors of pions and kaons, as well as those for Ki^ decays. We used 
the same framework to deal with masses, decay constants and form factors simulated on the lattice 
for 2-|-l dynamical fiavours (in the unitary limit where valence and sea quarks have the same masses). 
This allowed us to extract information on the pattern of three-flavour chiral symmetry breaking from 
lattice data on light-meson spectrum and form factors, with a signiflcant paramagnetic suppression 
of the quark condensate and the pseudoscalar decay constant in the three-flavour chiral limit. In the 
present article, we aim at confronting the possibility of significant paramagnetic suppression of the 
quark condensate with lattice data gathered on the topological susceptibility. In Sec. 3, we reconsider 
the derivation of the NLO expression for the topological susceptibility in the case of a suppression 
of the LO Nf = 3 chiral order parameters. We also discuss the implications of this formula for 
the determination of the quark condensate. In Sec. 4, we illustrate our approach by fitting results 



obtained by the RBC/UKQCD data on the topological susceptibility as well as the pseudoscalar 
meson spectrum. An important issue consists in the determination of the lattice spacings and the 
related discretisation errors, which we discuss in Sec. 5 before concluding. 

2 Resumined Chiral Perturbation Theory and 77 observables 

We have already discussed the salient features of the resummed xPT framework |47j . Let us briefly 
recall the basic ideas keeping in mind that no specific hypothesis is made upon the relative size of 
LO and NLO. 



• 



• 



First, a subset of observables is chosen for which a "good" overall convergence is assumed, i.e. 
the sum of LO and NLO terms is large compared to the remaining higher-order (HO) terms 
of the series (in relation with the correlators described by xPT): for example, the squared 
pion and kaon decay constants, F^, F"^ (related to the two-point axial correlator {A'^A'^)), the 
squared masses, F^M^, F^M^ {{d^A^^dyA^)) and the pion electromagnetic form factor, F-j^Fy 
{{A'^V^^A'^)). Their chiral expansion is performed in terms of xPT Lagrangian parameters up 
to NLO, keeping track of the HO which form the (small) remainders of the series. 

Furthermore, three fundamental LO quantities 

A-(3)^^^, .(3).^, ..=^. (8) 



FIMZ ' '"' F2 



m 



as well as HO remainders, are kept free. The first two quantities in eq. ([8]) are of particular 

relevance, since they express two main order parameters of Nj = 3 chiral symmetry breaking, 

the quark condensate and the pseudoscalar decay constant, in physical units. They also assess 

the saturation of the chiral expansion of F^M^ and F^ by their leading order. One of them 

can be traded for: 

_ X(3) _ 2m5o 

^^'^'~ Z(3) " M2 • ^^^ 

The ratio r measures the relative size of the quark masses in a framework where the strange 
quark is supposed to play a peculiar role in the chiral structure of QCD vacuum. 

• Finally, the parameters of the Lagrangian which appear at NLO, the so-called low-energy 
constants (LEG) Li, are expressed as functions of -^(3), Z{3), r, some conveniently chosen 
observables and the remainders. The relations thus obtained can be inserted into the chiral 
expansions of other observables, introducing additional remainders in an appropriate way. 

The important point is not to systematically trade leading-order terms for physical ones (for 
example 2niBQ is not traded for M^), as usually done, but to perform the replacement only when 
physical arguments indicate that the convergence of the series will be improved, and thus higher-order 
remainders will see their size reduced. For example, the nonanalytic structure imposed by unitarity 
should be located at the physical poles, thresholds. . . We will not further comment on this question 
here since we only will encounter tadpoles contributions in the following. A thorough discussion of 
an implementation of exact perturbative unitarity and exact renormalization scale independence and 
its illustration for nr] scattering can be found in ref. [38] q 

As indicated by eqs. ([3]) and ([5|), the propagation and coupling of the vr*^ and rj mesons will play 
a particular role for the analysis of the topological susceptibility. We will thus include a further 
element that was not discussed in detail in ref. [ITj, but presented in ref. [S], i.e., the decay constant 



■^Let us mention that there are sUght diflferences between refs. [48] and [47] concerning the treatment of the unitarity 
part, which affect the corresponding pieces in the form factors and scattering amplitudes, but not the observables 
considered here. 



and mass of the rj-meson which features as a propagating mode of the correlator x- In the physical 
case and in the isospin limit, we can write down relationships for the r] observables: 



F' 



J. ^ ±v±^ 



F^Z{3) + 8(r + 2)Y{3)M^ALi + i:i2r + 1)Y{3)M^AL5 



487r2 
2r + l 



(2r + 1) loe 



o 2 
M 



V 



o2 

Mk 



log 



o2 ■ 

Mk 

o 2 

M^ 



+ ^r} ^V ' 



16 



-^-F^M^XiS) + y (2r + l)(r + 2)y^(3)M;^AL6 



(10) 



(11) 



i-r4 



l)^y^(3)M^L7 + F^"M^^d^ 



Mp denotes the LO contribution to the mass of the pseudoscalar meson P = tt,K, rj. ALj = L^ — Li 
are scale- independent combinations of the NLO low-energy constants Li and chiral logarithms L[. 
Using NLO chiral series with expected good convergence properties, these combinations can be 
reexpressed in terms of pseudoscalar masses, decay constants, LO parameters and HO remainders, 
as discussed at length in ref. [4j^J and briefly recalled in App.|Xl The remainders (ir,,e^ are assumed 
to be small, with the estimate for the order of magnitude d^, e^ = 0{m'l) = 0(10%) [in practice, we 
will follow the dimensional estimates described in App. B in ref. |47)]. 

In ref. [47], we used the relationships between ^4^5 and F^, F^, as well as those between Lg^s and 
F!^M^,F'^M^ to express the NLO low-energy constants (LECs) in terms of measured quantities, 
HO remainders as well as the three LO parameters eq. ([8]). Following the same approach, we can 
invert eq. pip to reexpress the low-energy constant L7 in terms of masses and decay constants: 



[Yi3)fL, 



32(r - 1)2 M2 



3F^M^i - AF^Ml + F^Ml 



dco 



l?\e{r) 



+ d'] 



(12) 



where Frj is not measured, but can be computed using eq. (fTO]l . We have introduced the Gell-Mann- 
Okubo difference of masses and the equivalent for HO remainders: 



A 



GO 



3F^M^ - ^F^Ml + F^Ml 



F^M^ 



dco 



3F^M^ 



dr, 



7r 7r 



and the function of r: 



e(r) 



J2 



1' 



J FkMk X ^ 
ro = z — 1 ~ 



36, 



(14) 



so that e(r2) = 0, and e(ri) = 1 with n = 2{FkMk)/{F.„M.„) - 1 ~ 8. 

In the following, we will consider simulations with 2-1-1 dynamical flavours (m, m, fhs) and denote 
X the values obtained from lattice simulations (and X the corresponding value for physical quark 
masses). We introduce the ratios: 

T?) _ rh. 

(15) 



ms 
P= — 



m 
Q= -^ , 



in addition to the ratio of physical quark masses r and the chiral parameters arising in the LO 
Lagrangian in eq. ([8|). This yields the following expressions for the r] mass and decay constant for 



lattice simulations: 






F^Z{3) + 8pqr - + 2 Y{3)M^Ll{fi) + -pqr - + 1 Y{3)M^Ll{fi) 



r2rr/ 



~ 2 
^p,ry(3)logM^+F,%, 



327r 
pqr [2 
~3 

32 

1, 



+ 1 ) F^^MlXi?,) + ^^{pqrf f - + 1 V^ + 2 ) i^^(3)MiL^(/.) 



16/2 



+ "-^{pqrf ( - - ij Y\3)MtLr + y (^^ + 1 ) (p?r)2y2(3)M,^Li(/.) 



-#^->-(*5e- 



~ 2 ~ 2' 

, 1 \ Mk 4/2 \ , M„ 
2 - + 1 log ^ + - - + 1 loK ^ 



/^' 



9 Vg 



/^^ 



+ 



~ 2 



log^-iri + lVo.^-lf? + ll.or-" 



/^" 



3 \q 



ix^ 



9 V9 



~ 2 

o 

Af, 



+ F^M^d^ 



(16) 



(17) 



where L^ is given by eq. (|12|) . and equivalent expressions for -^4568 ^^"^ ^^ found in ref. [47] in terms 
of r, X(3), 2'(3), HO remainders, masses and decay constants (they are recalled in App. [X]). The HO 
remainders dr^, Crj are supposed to scale like m^ for simulations where the simulated strange quark is 
much heavier than the other ones. The LO contributions to the masses of the simulated pseudoscalar 
mesons read: 



~ 2 ~ 2 

M^ = pqrMlY{3) , Mk 



pqr 



1 



+ 1 MlY{3) , 



~ 2 



pqr 
~3~ 



+ 1 MlY{3) . (18) 



3 The topological susceptibility at next-to-leading order 



3.1 Diagrammatic analysis 

As discussed in the introduction, the topological susceptibility can be seen as the value at zero of 
the two-point winding number density correlator, or as the averaged winding number density. It 
can be computed in the effective theory by considering the whole correlator (as was done in ref. [13] 
to obtain the leading-order expression) or by considering only the effective potential for constant 
sources and determining its dependence on the vacuum angle 9 (like the derivation of the NLO 
expression in ref. [35]). In the latter case, we have performed the calculation keeping all orders in 
strong isospin-breaking (contrary to ref. [35j) and obtained in the xPT framework 
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(19) 
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Figure 1: Diagrams involved in the evaluation of the correlator x- The winding number densities uj are 
represented by the gray boxes, the propagating rj and tt'^ mesons by the solid lines attached to the sources, 
whereas the loops contain any of the eight Goldstone bosons. The black dot in diagram (e) corresponds to the 
usual four-point vertex, and the white square in diagram (f) to the NLO low-energy counterterms. 



where e is the LO vr'^ry mixing angle, defined as: 



tan2e 



%/3 



rrid 



m,, 



m.o 



m 



m = {niu + m^)/2. 



(20) 



This expression is indeed independent of the renormahsation scale /x and agrees with ref. [35] in 
the isospin limit, as expected. However, this result as it stands does not exploit the Ward identity 
eq. ([3]), which suggests strongly that one should identify the r/, ir^ poles in the chiral expansion of the 
correlator. 

It is straightforward to analyse the correlator x, eq. ([2]), diagrammatically in NLO xPT to identify 
the parts corresponding to propagation, as shown in Fig. [TJ The LO expression eq. ([5]) comes from 
the local terms in diagram (a) as well as the propagation part in diagram (b). At NLO, the local 
term gets contributions from diagram (c) , the coupling of the flavour-singlet mesons to the winding 
number density receives corrections from diagram (d) and the propagation of the mesons is corrected 
by diagrams (e) and (f), which are partly absorbed by the redefinition of the pole position, yielding: 



X{(1^ 



E 



Rp 



Ml 
■K",ri ^ 



+ S{q' 



(21) 



where Rp = |(0|a;|i-*)p is the relevant coupling to Goldstone bosons, and S is analytic up to the next 
thresholds (the first singularities being the cut from s > 9M^ and pole at M^,). For the topological 
susceptibility in the isospin limit and at next-to-leading order, one obtains the separation: 



X 



pole 






R^R, = mu;\rj)\' 



S = S{0) . 



(22) 



The vr'^ meson does not couple to the winding number density in this limit, due to G-parity. Up to 
next-to-leading order the pole residue and the analytical piece read: 
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(23) 
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+ -Bl [(2m + m,)2(L^(;x) + Ll,{ix)) + ilm" + m2)L^(^)] 



(24) 



It is easy to check that eq. (|22p coincides with eq. (|19|) in the isospin limit, up to higher-order terms 
in the chiral expansion. 

3.2 Resummed expression (no r/ pole) 

Like in the case of F!^M^, the three-flavour quark condensate arises at LO in eq. (|19p . whereas NLO 
involves the low-energy constant Lq, related to the violation of the Zweig rule in the scalar sector and 
enhanced by a factor of m^. This expression (or even its truncation at LO only) is used to determine 
the quark condensate, assuming the smallness of NLO corrections. In refs. [16H18|,H3Vt47j . we have 
argued that the pattern of Nf = 3 chiral symmetry breaking could be affected significantly by vacuum 
fluctuations of ss pairs, leading to the suppression of the quark condensate, the enhancement of Lq 
(and L4) and finally a numerical competition between LO and NLO contributions in Nj = 3 chiral 
expansions. Such a problem would occur for eq. ()19p . which explains why we want to analyse the 
topological susceptibility in the Re^PT framework. 

We consider now a lattice simulation with dynamical quarks of unphysical masses {rh,rh,rhs). If 
we do not isolate the contribution from the rj pole and simply reexpress eq. (|19p using eqs. ([8]) and 
P^ . we obtain: 
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I v'^° pole J 



(25) 



^yno pole ^ 



where d-^ denotes the remainder collecting higher-order contributions. If we replace -^6,7,8 by their 
expressions in terms of LO quantities and physical observables, eq. P^ and App. [XJ and if we 
consider a simulation at the physical point {q = l/r,p = 1), the topological susceptibility boils down 
to: 



X 



no pole 



F2M2 



2 2r + l 
1 



167r2' ' '' F2 



[y(3)]2 

M2 



1 - e(r) 



M2 



7r + 2 



9 



2r(2r + l) 2(r-l)2(r + 2) 



2r + l 



log 



^2 



+ 



(l+r)2 
r(2r + 1; 



^GO 



log 



Mk r + 2 



9r 



+ 



M" 



log- 



o 2 
^2 



L^(;u)(r + 2) + 3Lg(^) 



2r + l 



+ ... 



(26) 



where the ellipsis denotes higher-order remainders, and the combinations of chiral logarithms L\, 
which depend on Y[?>) and quark masses, are given in App. |X1 Eq. (j26p shows that in the physical 
case, the topological susceptibility has no sensitivity on the three- flavour condensate -^^(3). Indeed, 
up to 1/r-suppressed corrections (which are small since r > 15 [21 [3116]), eq. ([25|) combines the 
(LO) three-flavour condensate and the (NLO) low-energy constant Lg in the same way as F^M^, or 



equivalently the two- flavour condensate X{2): 



X{2) 



lim l^i^ 
m^o 2m 



F^M^X{3) + 16r[Yi3)fM^Ll{i,) + 



(27) 



which is determined from the study of vrvr scattering through Ki4^ decays [IHZj • As was discussed in 
ref. |16] . the two- flavour condensate can be significantly different from the three-flavour one, and its 
value is strongly correlated to the ratio of quark masses r through the three-flavour chiral expansion 
of F^M'^. This feature explains both the absence of X(3) and the r-dependence in eq. (f^Sj) . 

Therefore, the topological susceptibility close to the physical point {p = 1, r = 1/q large) exhibits 
a rather poor sensitivity to the three- flavour condensate ^(3). In a similar way to pseudoscalar 
masses and decay constants, simulations aiming at disentangling the pattern of three-flavour chiral 
symmetry breaking should be performed not only for quark masses with a hierarchy similar to the 
physical case, but also light (almost) degenerate masses with values between the physical mu.,md 
and mg. The (unphysical) region where simulations probe X(3) efficiently can be determined by 
expanding eq. ()25p in powers of 1/r: 



X 



no pole 



FIMI 



pqr 



q + 2 



X{?,)[l-p{2q+l)]+p{2q+l) 



3pq r2 



2{q + 2) 



+ 



(28) 



where the ellipsis denote chiral logarithms and HO remainders. It is clear from this expression that 
we need q = 0(1) rather than 0(l/r) to get a reasonable sensitivity to -^^(3): for instance, q = I (i.e., 
m = rhs yields [1 — 3p]X(3) +p[3 — (r2 — r)/(2r)], so that the coefficient of X(3) has a similar size 
to that of the constant term. This is sketched in a more quantitative way in Fig. [2] which compares 
the topological susceptibility at X{?>) = and X{?>) = 1 in some illustrative cases Q The comparison 
between the two values indicates the sensitivity to the size of the three-flavour quark condensate, 
which tends to decrease the value of the topological susceptibility. We see that the contributions 
from X{3) are at least twice as small as the remaining contributions, and they in particular might 
be hidden in the uncertainties if q is small (e.g., of 0(l/r) if the simulated quark masses are tuned 
to be close to the physical ones). 

Another way of escaping the poor sensitivity of the topological susceptibility to the three-flavour 
quark condensate consists in simulations involving significant isospin breaking, which can be easily 
implemented in particular for twisted-mass fermion actions |32y5U) . An expression for the topological 
susceptibility on the lattice similar eq. ()25p can be written, involving the two quantities qu = rhu/rhs 
and qd = rhd/rhs instead of q: 



X 



no pole 



F^M^ pquqdT 
2 qu + qd + quQd 

pr{l + qu + qd) 



(29) 



X(3) 1 



r + 2 



e{r)pr 



^uQd 



2(gn + qd + QuC, 



+ 



l + qu + Qd 

r + 2 



pr{l + qu + qd) 

r + 2 



QpquQdr 



2{r - lY{qu + qd + quqd) 



A, 



GO 



+ 



where we expressed the NLO LECs ^6,7,8 using App.|A](we have checked that including the md — rriu 
difference in the expressions of App. |X]would not modify our conclusions). The ellipsis denotes chiral 
logarithms and HO remainders. The dependence of x on the two quantities qu and qd is illustrated 
in Fig. El As expected, x{qu,qd) = x{qd-,qu) vanishes for g„ = or g^ = since the theory becomes 
independent of 6 when one of the quark masses vanishes. The value of the topological susceptibility 



^For this comparison, we assume that the physical value of F^ is known in order to compute the Gell-Mann-Okubo 
contribution Ago to L-j. In the following sections, we will not make this assumption and we will use eq. (|10|) to 
determine _F„. 
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Figure 2: The topological susceptibility x^° p° as a function of the ratio of simulated quark masses q for 
X(3) = (solid lines) and X{3) = 1 (dashed lines) for two different values of p (related to the simulated 
strange quark mass). For illustration, the remainders are set to be zero, F(3) = 0.8, r = 26, Fk/Fj, = 1.19 



andF.,^ ^130 MeV. 



increases when one of the two quark masses increases, the other one being kept fixed. As clearly 
seen in Fig. [3l the contribution to the topological susceptibility independent of X{'i) increases faster 
than the one proportional to ^(3), so that a scan through values of {mu,rh(i) at fixed rhg could 
help to determine unambiguously the contribution (and thus the value) of the three-fiavour quark 
condensate. 



3.3 Resummed expression (identifying the r] pole) 

Identifying the poles corresponding to the propagation of the pseudoscalar mesons is of particular 
relevance within Resummed xPT. In this framework, smaller HO remainders are expected when 
singularities of the expansion (and in particular poles) are correctly located at the physical positions. 
Let us notice that a similar discussion takes place when one considers the two-point correlators of 
the axial current and/or pseudoscalar density, where the propagation of a single Goldstone boson 
occurs also at LO. Indeed, in refs. p^Hl8lH3tii7] . we have considered the decay constants and masses 
as they could be obtained from the NLO expression of (A/^iA^) where the propagators are explicitly 
expressed in terms of the physical masses, considering either the correlator at (for F^) or its residue 
at the pseudoscalar pole (for F^M^). 

Thanks to our diagrammatical analysis of the two-point correlator x iii xPT, we can easily identify 
the contribution from the propagation of the r] meson in eq. (|25|) . and we obtain the expressions for 
the pole residue and the analytical piece arising in eq. (|22p : 



R 



{pr{q - l)f 



27 



F^M:^Xi3)Y{3) + -ipr)\l 



''M^Y{3y 



(2q + l)Ll + ^Ll 



32 



-^{prrY{3rM^[{Ll + Lr){q - l)^(2g + 1) + Ll{q' - l){q - 1) 



^(H^Jg; 1)^6^(3)3 



8647r2 



~ 2 



Qq' log ^ - (1 + qf log ^ + ^(g - 4)(, + 2) log ^ 



(30) 
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Figure 3: The topological susceptibility x"° ^°^° (in units of GeV^) as a function of g„ — rhu/fris and qd = 
rhd/iTis for X(3) — (solid) and X^ = 1 (dashed with colour levels). For illustration, the remainders are set 
to be zero, p = 1, Y{3) = 0.8, r = 26, Fk/F^ = 1.19 and F,-, ^130 MeV. 



P-!^±Af^^MIX{?,) + l{prfY{2,fMtm + l)\Ll + L,) + {2q^ + l)Lg] (31) 

(pr)^ 
~ 5767r2 



2ji/r4 



Y{?,YM'* 






2 +-(2 + g)^log 



~ 2 



/^^ 



where the counterterms can be expressed in terms of r, X(3), .^(3) and the remainders. The topolog- 
ical susceptibility then reads: 

^ +5 + F°%Poic. (32) 



^pole 



M2 



At this stage, one has still to discuss the direct remainder xdx which corresponds to higher-order 
terms, and has been added in eq. ()32p as well as in eq. ()25p . One wants d^ to start at NNLO, and 
thus to be 0{m'i). This expectation can be checked by considering various chiral limits sending m 
and/or fhs to zero. As already explained previously, the topological susceptibility x should vanish 
in the limit where one of the quark masses goes to zero. In our framework, this can be translated as 
X — 7" for the three chiral limits: 



m.o 



m 



■> 0, ?fi fixed: p ^- 0, pq fixed, 
0, rhs fixed: q ^ 0, p fixed. 



• m, rhs — ^ 0, rh/rhs fixed: p — t- 0, (7 fixed. 

By inspection, one can check that the sum of LO and NLO in eq. ()25p does vanish in these three 
limits, whereas eqs. ()30p -()3ip do vanish in these limits only up to nonzero higher-order corrections. 
Therefore, if we consider the expression for the topological susceptibility not singling out the ry-pole, 
eq. (p5]) . we can add a higher-order remainder of the form x^" P°^''(i^no poie with d^ \>oe _ q(j^t^ thai 
has no singularities when any fhq — )• 0. In the situations that we consider (one of the masses much 
larger than the other ones), we expect d^ ^°'^ to be dominated by ml contributions. 
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Figure 4: The topological susceptibility ;^p°'° as a function of the light quark mass. The black [red] points 
correspond to the (24)"^ x 64 volume [{32)^ x 6AJ, using the estimates of lattice spacings quoted in ref. \S(Jf . The 
dashed curves correspond to the best fit using ;^p°'° with finite-volume effects, as indicated in the last column of 
table 1 (removing finite- volume effects would lead to almost identical curves). We have not included the lightest 
point in our analysis, as discussed in the text. As a reference, we have also indicated the same variation when 
the strange quark mass is set to its physical value. 



On the other hand, if we single out the ry-pole contributions following eq. (|25|) . we would need 
d^ ^ to become singular in the chiral limits considered before, so that x^'^^'^d^poic does not vanish and 
cancels the non- vanishing value of eqs. ()30p -()3ip coming from higher-order terms. A more satisfying 
solution consists in subtracting these higher-order pieces from the HO remainder, so that d^ does 
not exhibit singularities in any of the chiral limits. In other words, using eqs. (|30| ) - (|31| ). we write the 
topological susceptibility as: 






X 



pole 



C{m,ms 



C(0, rhs) - C{m, 0) + xP°''dvPoic . 



(33) 



One can check that the subtracted terms C{0,rhs) = 0{rh'^) and C{m,0) = 0{fh?) are indeed HO 
terms. With the definition eq. ()33p . which will be used in the following, d^°^ is regular in the chiral 
limits described above, and can be considered as 0{m1) for the simulations considered here. 

Finally, one should notice that R and S involve again not only the three-flavour condensate X(3) 
but also NLO LECs, and in particular Lg. This dependence is equivalent to that of eq. ([25]) up 
to higher orders in the expansion of the quark masses. It is therefore of no surprise that we find 
only small numerical differences between eq. (|25p and eqs. (J30p -(j3i p in the range of quark masses 
of physical interest, and that the previous conclusion concerning the potential of the topological 
susceptibility to determine the three-flavour quark condensate still applies in this setting. 



3.4 Finite-volume effects 

Two different lattice artefacts can affect the results of simulations before considering the continuum 
limit. The first corresponds to the finite size of the lattice spacing, which can be included in principle 
in the chiral expansion through spurious that depends on the implementation of the fermion action 
(see ref. |51] and references therein). For the moment, we will not include this effect as we will 
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consider data with good chiral properties and only 0(0^) artefacts, even though this issue is naturahy 
of interest [26] . 

A second effect comes from the finiteness of the volume used for lattice simulations, inducing 
finite- volume modifications of the chiral expansions. As discussed in ref. |19 y52f[55] . the NLO finite- 
volume effects for chiral perturbation theory amount to a modification of the chiral logarithms: 



Mp Mp 



J_ V— - — 



2L3 ^ ojp 



(34) 



where i- G 2it/L x Z^ and wp = \l 1"^+ Mp- The summation over (Fourier conjugates of) the three 
spatial directions comes from the quantization of momenta due to the periodic boundary conditions 
on the lattice box (making it a torus in practice). We consider simulations where the time component 
is significantly larger than the spatial components, and consider only the finite- volume effects related 
to the latter |53H55j . One can make contact with the infinite- volume logarithm through the function: 



ap_ 
L3 



o2 

Mp 
87r2 



o 2 

log -^ + r.p , 



(35) 



o2 



where Hp = ^i/2{L,M p) was introduced in ref. [53] . An alternative definition of Hp was proposed 
in ref. [l5] and is rediscussed in App. O 

We will take into account these corrections in the following analysis for all the quantities of 
interest. For masses and decay constants, the corrections can be read directly from the expressions 
in ref. |47j . For the topological susceptibility, we have the following correction from finite- volume 
effects to the expression without singling out the r/-pole: 



X 



no polc^^^ ^ ^no pole^^^ 



1 



n3)M|j^ 



and when one singles out the r/-pole: 
R{L) 



i?(oo) + M%^y(3)^M.^ 



S{L) 



5(oo) 



{pr) 



54 



y(3)Mj 



6H^ + 4((? + l)Hx + -(2g + l)H^ 



3gH^ - (1 + qy^K + 3(9 - 4)^'? 



6(7H^ + 4(g + l)Hi^ + -(2 + (7)H^ 



(36) 

(37) 
(38) 



We neglect any dependence of the higher-order remainders on finite-volume effects, effectively iden- 
tifying the HO remainders at finite and infinite volumes. 

4 A first series of fits to lattice data 

We want now to exploit the topological susceptibility to improve our determination of the pattern 
of three-fiavour chiral symmetry breaking. We have however seen that the usual setting of lattice 
simulations is not appropriate, as the topological susceptibility is then driven by the two-flavour 
condensate. If we want to extract information on the three-flavour pattern, we need either to make 
simulations away from the physical case (as discussed in Sec. 13. 2p . or to supplement the topological 
susceptibility with other sources of information. We will now take this second option, which is also 
required due to the number of parameters involved in the NLO expression of the topological sus- 
ceptibility (in particular the parameters of the LO chiral Lagrangian and HO remainders). We will 
therefore include pseudoscalar masses and decay constants in our considerations and test the com- 
patibility of the topological susceptibility with the pattern of three-flavour chiral symmetry breaking 
obtained from the latter data. 
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There are several lattice calculations of the topological susceptibility in the literature, e.g. refs. 
130] . As an illustration of our approach, we will focus on the ones from RBC/UKQCD following our 
work, ref. [17], as this collaboration provides all the details (masses, decay constants and topological 
susceptibility) necessary for our analysis. Since the publication of ref. |47] . new data from this 
collaboration have been issued with a new volume and different quark masses, given in ref. [30] . 
and we will use them as a reference in the coming sections. We follow the approach of ref. [37] and 
perform a fit to the two sets of data in the 24^ x 64 (a~^ = 1.78(3) GeV) and 32^ x 64 volumes 
(a^^ = 2.28(3) GeV) with and without including finite volume effects. In this section, the lattice 
spacings are fixed to the central values quoted in ref. |3D] without uncertainty. 

The fits include the data collected in App. |Bl i.e. pion and kaon masses and decay constants as 
well as the topological susceptibility (unfortunately, no lattice data on the t] meson is available, even 
though it would improve the analysis of the topological susceptibility in a significant way) . We build 
a x^ depending on the following parameters: 

• the three leading-order parameters r,X{3),Z{3), 

• a reference ratio between a simulated strange quark (chosen conventionally to be for the 24^ 
simulations) and the physical strange quark mass, 

• the HO remainders associated with the pion and kaon masses and decay constants (denoted 
d, e, d' , e'), 

• the ratio Fk/F-,^ (on the other hand, we set F^^ = 92.2 MeV), 

• if the topological susceptibility is included, the corresponding HO remainders for the rj mass 
and decay constant (d,,,e^) as well as the one for the topological susceptibility (d^,). 

HO remainders are restricted to remain within a range based on a naive dimensional analysis, as 
described in ref. [47J. At this stage, we include no discretisation error effects in our three-flavour 
chiral expansions, but we will come back to this issue in Sec. [S] 

As shown in Tab. [U the pattern of A'^^ = 3 chiral symmetry breaking with low quark condensate 
and decay constant observed in that reference is confirmed by this new analysis, whether we include 
or not the topological susceptibility in our fit. The outcome of the fit is thus mainly driven by the 
spectrum of pseudoscalar mesons, but the quality of the agreement is not modified by the inclusion 
of the topological susceptibility, which thus exhibits a good compatibility with the rest of the fit and 
is consistent with the pattern of chiral symmetry breaking described in ref. [47] . As in our previous 
work, L4 and Lq do not show any sign of Zweig suppression and the competition between LO and 
NLO in three-flavour chiral expansions is clearly seen. We obtain values for the Nf = 2 chiral order 
parameters in agreement with expectations from two-flavour chiral perturbation theory |13j as well 
as experimental information on vrvr scattering, such as that from Ki4^ decays [IH7]. The values of 
•£3 and £4 given there can be compared to the one obtained from the RBC/UKQCD collaboration, 
ref. [30]: 4 = 2.82(16), £4 = 3.76(9) in infinite-volume xPT and 4 = 2.57(18), £4 = 3.83(9) in finite- 
volume xPT. One can also recall the value quoted by the Flavour Lattice Averaging Group |15] for 
£3 = 3.2(8) (no value was quoted for ^4 in this reference). Finally, one notices that singling out 
the r/-pole or not in the expression of the topological susceptibility does not modify significantly the 
analysis - the values obtained with eqs. ()25p and ()30p - ()3ip are very close numerically for the ranges 
of parameters considered here. 

Focusing on the results including x^ with finite-volume effects, we obtain the following conver- 
gence at Xmin with the relative LO, NLO and HO contributions: 



^pole 

A, 



0.65 + 0.40 - 0.05 , F^M^ 

0.47 + 0.56 - 0.03 , F^^Ml 

0.40 + 0.57 + 0.03 , F^M^ 
0.22 -F 0.63 + 0.15, 
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0.35 + 0.79-0.14, 
0.24 + 0.86-0.10, 
0.23 + 0.80-0.03, ^ ^ 



Fit name 


Al 


A2 


A3 


AA 


A5 


X 


No 


No 


^nopolc 


^pole 


^pole 


Finite volume 


No 


Yes 


No 


No 


Yes 


r 


23.0 ±0.7 


23.4 ±0.7 


23.0 ±0.7 


23.0 ±0.7 


23.4 ±0.7 


X(3) 


0.33 ±0.10 


0.34 ±0.05 


0.33 ±0.06 


0.33 ±0.06 


0.34 ±0.05 


y(3) 


0.49 ±0.15 


0.53 ±0.07 


0.49 ± 0.08 


0.49 ± 0.08 


0.53 ±0.07 


Z(3) 


0.68 ±0.03 


0.65 ±0.03 


0.68 ± 0.03 


0.68 ± 0.03 


0.65 ±0.03 


Fk/F^ 


1.17 ±0.01 


1.17 ±0.01 


1.17 ±0.01 


1.17 ±0.01 


1.17 ±0.01 


Rem. at limit 


d 


d 


d, e, d^ 


d, e, d^ 


d,e.^,d^ 


ms{2A^)/ms 


1.12 ±0.02 


1.12 ±0.02 


1.12 ±0.02 


1.12 ±0.02 


1.12 ±0.02 


ms{2 GeV)[MeV] 


98.1 ±1.7 


98.4 ±1.8 


98.0 ±1.7 


98.0 ±1.7 


98.4 ±1.8 


m(2 GeV)[MeV] 


4.3 ±0.1 


4.2 ±0.1 


4.3 ±0.1 


4.3 ±0.1 


4.2 ±0.1 


t}J^{2 GeV)[MeV] 


186 ± 19 


211 ±7 


186 ±11 


186 ±12 


189 ± 10 


So (2 GeV)[GeV] 


1.12 ±0.34 


1.22 ±0.17 


1.11 ±0.21 


1.12 ±0.21 


1.22 ±0.17 


Fo[MeV] 


75.9 ±1.3 


74.4 ±1.4 


75.9 ±1.3 


75.9 ±1.3 


74.5 ±1.4 


F^[MeV] 


— 


— 


117 ±16 


128 ± 10 


124 ±6 


X-10^[GeV4] 


— 


— 


0.51 ±0.01 


0.49 ± 0.02 


0.50 ± 0.02 


m,(323)/m. 


1.08 ±0.02 


1.09 ±0.02 


1.07 ±0.02 


1.08 ±0.02 


1.07 ±0.02 


U{^^) ■ 10^ 


1.12 ±0.46 


0.65 ± 0.45 


1.13 ±0.30 


1.12 ±0.30 


0.60 ±0.43 


i5(^) • 103 


2.13 ±0.78 


2.05 ± 0.40 


2.14 ±0.53 


2.13 ±0.53 


2.04 ±0.39 


L(,{^i) ■ 10^ 


3.00 ±3.05 


2.55 ±0.89 


3.13 ±1.44 


3.10 ±1.42 


2.52 ± 0.88 


L7(^) • 103 


— 


— 


-2.60 ±1.02 


-1.82 ±0.47 


-1.70 ±0.40 


L8(/i) • 103 


4.12 ±2.74 


3.39 ±1.13 


4.17 ±1.79 


4.12 ±1.77 


3.35 ±1.11 


X{2) 


0.90 ±0.01 


0.90 ±0.01 


0.90 ±0.01 


0.90 ±0.01 


0.90 ± 0.01 


Y{2) 


0.99 ±0.01 


1.00 ±0.01 


0.99 ±0.01 


0.99 ±0.01 


1.00 ±0.01 


Z{2) 


0.91 ±0.01 


0.91 ±0.01 


0.91 ±0.01 


0.91 ±0.01 


0.90 ± 0.01 


T}/^{2 GeV)[MeV] 


260 ±3 


261 ±3 


260 ±2 


260 ±3 


261 ±2 


B{2 GeV)[GeV] 


2.26 ±0.07 


2.31 ±0.06 


2.26 ±0.07 


2.26 ±0.07 


2.31 ±0.06 


F[MeV] 


88.0 ±0.4 


87.8 ±0.2 


88.1 ±0.2 


88.0 ±0.22 


87.8 ±0.14 


h 


-1.6 ±1.4 


-0.5 ±0.9 


-1.6 ±1.2 


-1.6±1.2 


-0.5 ±0.9 


U 


3.0 ±0.3 


3.2 ±0.1 


3.0 ±0.1 


3.0 ±0.2 


3.2 ±0.1 


Y^IY.^ = S(2)/S(3) 


2.72 ±0.78 


2.64 ±0.35 


2.74 ±0.46 


2.72 ± 0.46 


2.62 ±0.35 


B/Bo = B{2)/B{3) 


2.02 ±0.56 


1.89 ±0.23 


2.03 ±0.32 


2.02 ±0.33 


1.89 ±0.22 


F/Fo = F(2)/F(3) 


1.16 ±0.02 


1.18 ±0.02 


1.16 ±0.02 


1.16 ±0.02 


1.18 ±0.02 


xVn 


28.7/11 


18.7/11 


29.2/12 


29.4/12 


19.2/12 


Gaussian equiv. 


3.0(7 


l.So- 


2.9a 


2.9a 


1.7a 



Table 1: Results of fits performed on the data from RBC/UKQCD collaboration on pseudoscalar masses and 
decay constants, including or not the topological susceptibility '30]. In two cases, finite-volume effects are taken 
into account. In all cases, we considered only data with light pions and only statistical errors are shown. The 
LECs are given at the scale n = Mp. 

showing that the global convergence is good (small HO remainders), but the LO and NLO contribu- 
tions are indeed competing numerically, confirming the results obtained in ref. |47) . 

However, fitting the RBC/UKQCD data in the two volumes, we obtain a rather poor fit of 
X^/d.o.f.=29.4/12 without finite-volume effects, which can be seen as a 2.9a discrepancy in a naive 
(purely Gaussian) statistical interpretation. The fit improves when one includes finite- volume effects, 
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Figure 5: The pion decay constant as a function of the light-quark mass. Blue [red] points correspond to the 
data for the (24)'^ x 64 volume /(32)'^ x 64/, using the estimates of lattice spacing given in ref iSO^ . Purple 
[orange] points correspond to the same points, using our own estimates of lattice spacing discussed in Sec. ^1 
The dashed curves indicate our best fit in each case (015 and Bb), as given in the last columns of tables 1 
and 2 respectively. They are interrupted for A'l-^L < 2 where our description of finite-volume effects becomes 
unreliable. The solid curves indicate the corresponding variations for the physical value of the strange quark 
mass, without finite-volume effects. In both cases (estimate of the lattice spacings from ref. iSOjj or Sec. \^, the 
position of the physical point (jTjphys, -FV) *s indicated with a square. 



getting down to a 1.7a effect, as can be seen by comparing the fits Al and A2 (or AA and A5) in 
Tab. [TJ Indeed, even thougli these effects are rather smah in the ranges of quark masses considered 
here, they tend to bring the various quantities in better agreement with the lattice data (more on 
this issue can be found in App.[C]). The main contribution to Xmin comes from F^ (we wih come back 
to this issue in the next section) whereas the topological susceptibility contributes only marginally. 

In Fig. m we illustrate the results of the fit for the topological susceptibility as a function of the 
simulated light-quark mass. Finite-volume effects have little impact, since at large rh, the product 
of Mt^L is much larger than 1 and thus finite-volume effects are small, whereas at small rh, the 
topological susceptibility (with or without finite volume effects) goes to linearly. One notices also 
the linearity of the three curves, related to the fact that q remains small (< 0.4) in this range of m, so 
that X in eq ()25p or ()30p -()3i p can be expanded in powers of q with only small 0{q^) corrections. As 
is clear from eq. (|25p , the slopes of these curves at zero are not directly related to the three- flavour 
quark condensate, as it involves also the NLO LEG Lg as well as chiral logarithms (once again, this 
feature is more related to the two- flavour quark condensate). 

The authors of ref. [30] noticed that their value for the topological susceptibility obtained from 
the simulation with the lightest quark mass may suffer from "long auto- correlations in Q that are 
not well resolved in [their] finite Markov chain of configurations". Indeed, the uncertainty attached 
to this particular point is much smaller than the other ones, which might indicate underestimated 
systematic effects. In addition, one notices that it is much lower than the other points, so that it 
would be very difficult to draw a smooth curve going through all the points as well as the origin (as 
required by the vanishing of the topological susceptibility when one quark mass is set to zero). For 
these reasons, we did not include this data point in our analyses. We have however checked that 
including this additional point leads only to a worse \^ without affecting the outcome for the LO 
parameters of the chiral Lagrangian in a significant way. 
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5 The role of lattice spacing 

Our approach, allowing for a numerical competition between LO and NLO contributions to three- 
flavour chiral series, has provided a good, but not completely satisfying, fit of masses, decay constants 
and topological susceptibility from the RBC/UKQCD data. This is illustrated in the lower part of 
Fig. [5l where we plot F-n- as a function of rh (as given by the equivalent of eqs. ()16p -()17 p for F^ 
and F^M^, see Sec. 5.1 in ref. [47J). The (dashed) curves correspond to our best fit (last column in 
Tab. [T]) and the solid line indicates the dependence of F-n- for a physical ttt-s in an infinite volume. 
Indeed, in spite of this broad agreement, we notice that we get a large contribution to Xmin from the 
Ft^ data points. This is a reminder of the problem encountered in ref. |30j, where neither a chiral 
nor an analytic extrapolation formula was able to accommodate the observed dependence of F.,^ on 
fh with the physical point (w,phys,F^). Our formalism can include both pieces of information with 
a reasonable Xmin' ^^^ ^^ '^^^ improve the latter by letting the physical value of F^^ vary as a free 
parameter. We obtain the results indicated as fit Bl in the first column of Tab. [2] with a very low 
value of -Ftt = 86.4 MeV, in agreement with the results in ref. |30| . 

This discrepancy between lattice and physical values of F.,^ hints at a more general issue concerning 
the determination of the absolute scale of lattice quantities, which is obtained by converting the 
lattice spacing into physical units. In ref. [30j, it was determined by identifying "scaling trajectories" 
corresponding to lines in the {a,rh,ms) space with fixed Mj^/Mq and Mk/Mq ratios. An iterative 
interpolation method was used to reach values of Mt^,Mk,M^ corresponding to physical values of 
their ratios, which was identified as the physical point for the quark masses. The lattice spacings were 
then determined by requiring that 1/a = 1. 672/ (aM^) GeV where 1.672 GeV is the physical mass 
of this baryon and oMq is the mass of the as measured in lattice units. The actual interpolation 
was performed through two kinds of interpolating formulae for the hadron masses in terms of quark 
masses, either based on NLO three-flavour xPT or an analytic (polynomial) ansatz, fitted to partially- 
quenched data (where sea- and valence-quark masses are different). This led to values of the lattice 
spacing that were compatible and quoted with an accuracy at the level of a few percents. 

However, such determinations based on the quark-mass dependence of light-meson quantities 
might be affected significantly if one takes into account a numerical competition between LO and 
NLO contributions to three-flavour chiral series. In this particular case, one should consider at the 
same time the chiral expansions of FpMp and Fp and include HO remainders, to determine the 
dependence of M,r and Mk on the quark masses. As far as the mass is concerned, we will follow 
ref. [30] and take a linear dependence on the quark masses |j: 

Mn = Mq + ci{ms - rus) + C2{rh - m) , (40) 

inspired by the analysis of RBC/UKQCD. A fully consistent treatment should include a treatment 
of the baryon masses in our RexPT framework, extending eq. (j40|) to chiral logarithms and HO 
remainders. Such an analysis, beyond the scope of the present article, is under way [2|)j (in the 
present case, we have checked that adding a quartic term in eq. (pUj) does not change the results 
discussed below). 

In addition, since we are interested in effects related to lattice spacing, we should also consider 
discretisation errors, which could reach 10-15% in the analysis of ref. [3^. We follow the latter 



^The n mass was analysed as a function of mu,d in Nf = 2 chiral perturbation theory in ref. [56]. The Q,~ field is an 
isoscalar under SU{2), which prevents it from interacting with other baryons (contrary to other hyperons like E, E, A) 
in two-flavour xPT. This makes the dependence of AIn on M,r much simpler than for other nucleons, with a constant 
term supplemented with quadratic and quartic terms in Af^ (the quartic term including a chiral logarithm). 
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analysis and include only leading-order discretisation effects affecting the decay constants |j: 

Fl{a,V,m,ms) = F^{0,V,m,ms) + F^Z{3)cFy , (41) 

F|(a,y,m,m,) = FJ^{0,V,m,ms) + F^Z{3)cFj,a\ (42) 

where the correction term is defined with respect to the leading-order term in the chiral expansion. 
We are not in a position to perform the same joint determination of lattice spacings and quark 
masses as the RBC/UKQCD collaboration to include the dependence of the pion and kaon masses 
on the light-quark mass inferred from RexPT from scratch. However, we can perform a combined fit 
of the pion and kaon masses and decay constants, as well as the Q mass (collected in App. |B]), with 
the following parameters: 

• the three leading-order parameters r,X{3),Z{3), 

• a reference ratio between a simulated strange quark (chosen conventionally to be for the 24^ 
simulations) and the physical strange quark mass, 

• the HO remainders associated with the pion and kaon masses and decay constants (denoted 
d, e,d' , e'), 

• the ratio Fk/F.,^ (on the other hand, we set -F^ = 92.2 MeV), 

• if the lattice spacings are left free, the two effective constants ci, C2 for the fi mass, 

• if discretisation errors are included, the two effective constants cp^, cpj^ for the decay constants, 

• if the topological susceptibility is included, the corresponding HO remainders for the r/ mass 
and decay constant (d,,,e^) as well as the one for the topological susceptibility [d^). 

Our approach is not very different in its spirit from the "combined scaling and chiral fitting" per- 
formed in ref. [3U], up to the following modifications: we include information on the values of the 
masses and decay constants at the physical point, we consider FpMp and Fp rather than Mp and 
Fp, we use three-flavour Re^PT rather than two-flavour expansions to perform the interpolation 
of the data, we include the presence of HO remainders, we do not include partially-quenched data, 
we flx at the same time m, rus and the lattice spacings in contrast with the two-step procedure in 
ref. [30] (physical masses flrst, lattice spacing afterwards). 



^In principle, one should consider all the terms coming from discretisation effects and due to the breaking of chiral 
symmetry, and also add correction terms for the masses. This would however increase the number of free parameters 
and lead to fits with a poorly stability, due to flat directions in the subspace of correction terms and the limited number 
of data points. 
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Fit name 


Bl 


B2 


B3 


BA 


B5 


X 


No 


No 


No 


No 


^polc 


Finite volume 


Yes 


Yes 


Yes 


Yes 


Yes 


Lattice spacing 


Fixed 


From Mq 


Fixed 


From Mq 


From Mq 


0{a'^) corr 


No 


No 


Yes 


Yes 


Yes 


F.[MeV] 


86.5 ±2.2 


92.2* 


92.2* 


92.2* 


92.2* 


r 


25.6 ±1.1 


26.1 ±1.1 


23.3 ±0.8 


25.6 ±1.1 


25.6 ±1.0 


X(3) 


0.40 ±0.04 


0.54 ±0.05 


0.35 ±0.07 


0.55 ±0.04 


0.55 ± 0.04 


F(3) 


0.79 ±0.10 


0.86 ± 0.08 


0.52 ± 0.05 


0.91 ± 0.06 


0.91 ±0.07 


Z(3) 


0.50 ±0.07 


0.62 ±0.05 


0.68 ± 0.07 


0.60 ± 0.05 


0.60 ± 0.04 


Fk/F^ 


1.23 ±0.03 


1.25 ±0.03 


1.17 ±0.02 


1.25 ±0.02 


1.25 ±0.02 


Rem. at limit 


— 


d 


e 


d,e 


u, e, Ufq, e^q, (t^ 


m,(243)/m. 


1.14 ±0.03 


1.62 ±0.16 


1.12 ±0.03 


1.71 ±0.13 


1.71 ±0.13 


a-i(243)[GeV] 


1.73* 


2.14 ±0.12 


1.73* 


2.02 ±0.08 


2.19 ±0.08 


a-i(323)[GeV] 


2.28* 


2.76 ±0.13 


2.28* 


2.82 ± 0.09 


2.81 ±0.09 


ci 


— 


9.08 ± 0.44 


— 


8.98 ±0.37 


8.98 ±0.37 


C2 


— 


1.42 ±0.43 


— 


1.44 ±0.44 


1.44 ±0.43 


CF. [GeV2] 


— 


— 


0.13 ±0.14 


-1.09 ±1.19 


-1.09 ±1.19 


CF^ [GeV2] 


- 


- 


0.13 ±0.17 


0.09 ± 0.48 


0.09 ± 0.48 


m,(2 GeV) MeV 


96.5 ±2.2 


83.9 ±8.1 


97.9 ±2.2 


81.2 ±6.0 


81.2 ±6.0 


m(2 GeV) [MeV] 


3.8 ±0.2 


3.2 ±0.4 


4.2 ±0.1 


3.2 ±0.3 


3.2 ±0.3 


So/^(2 MeV) [GeV] 


197 ±7 


240 ± 16 


191 ±2 


243 ± 12 


243 ± 12 


5o(2 GeV) [GeV] 


2.03 ±0.32 


2.60 ±0.52 


1.21 ±0.12 


2.81 ±0.39 


2.81 ±0.38 


Fo[MeV] 


61.5 ±5.2 


72.8 ±2.8 


76.0 ± 3.7 


71.3 ±2.4 


71.3 ±2.4 


i^,[MeV] 


— 


— 


— 


— 


123 ±3 


^polc.iQ4[QgY4] 


— 


0.33 ±0.03 


0.41 ± 0.03 


0.33 ± 0.03 


0.43 ± 0.03 


m,(323)/m. 


1.09 ±0.03 


1.52 ±0.15 


1.08 ±0.03 


1.61 ±0.12 


1.61 ±0.12 


U{^i) ■ 10^ 


0.28 ±0.29 


-0.11 ±0.14 


0.99 ±0.17 


-0.08 ± 0.09 


-0.08 ± 0.09 


L5(/i) • 10^ 


1.62 ±0.19 


2.13 ±0.21 


2.03 ±0.43 


2.00 ± 0.20 


2.01 ±0.20 


L6(/x) • 103 


0.47 ±0.35 


0.09 ±0.10 


2.53 ±0.79 


0.05 ± 0.07 


0.04 ± 0.07 


L7(/") • 103 


— 


— 


— 


— 


-0.25 ±0.14 


Lsifi) ■ 103 


1.08 ±0.38 


1.13 ±0.20 


3.42 ±0.88 


1.04 ±0.14 


1.04 ±0.09 


X{2) 


0.89 ±0.02 


0.89 ±0.01 


0.90 ±0.01 


0.88 ±0.01 


0.88 ±0.01 


Y{2) 


1.03 ±0.02 


1.03 ±0.01 


0.99 ±0.01 


1.04 ±0.01 


1.04 ±0.01 


Z{2) 


0.86 ±0.02 


0.86 ±0.01 


0.91 ±0.01 


0.86 ±0.01 


0.85 ±0.01 


Si/3(2 GeV) [MeV] 


258 ±4 


248 ±2 


261 ±3 


285 ±7 


285 ±8 


B{2 GeV) [GeV] 


2.68 ±0.16 


3.13 ±0.42 


1.90 ±0.18 


3.18 ±0.28 


3.18 ±0.27 


F[MeV] 


80.2 ±2.9 


85.5 ±0.50 


88.0 ±0.4 


85.2 ±0.4 


85.2 ±0.4 


4 


4.0 ±1.5 


4.5 ±0.9 


-1.0 ±0.9 


4.6 ±0.7 


4.6 ±0.7 


^4 


4.2 ±0.4 


4.7 ±0.3 


3.1 ±0.3 


4.9 ±0.3 


4.9 ±0.3 


S/Eo = S(2)/S(3) 


2.24 ±0.22 


1.66 ±0.14 


2.54 ±0.02 


1.61 ±0.11 


1.61 ±0.11 


S/Bo = B{2)/B{'i) 


1.32 ±0.17 


1.20 ±0.10 


1.90 ±0.18 


1.13 ±0.09 


1.13 ±0.08 


F/Fo = F(2)/F(3) 


1.30 ±0.08 


1.17 ±0.05 


1.15 ±0.06 


1.19 ±0.04 


1.19 ±0.04 


xViv 


12.1/10 


17.3/12 


18.1/9 


16.0/10 


16.2/11 


Gaussian equiv. 


l.lcr 


1.5cr 


2.1cr 


1.6cj 


1.5(T 



Table 2: Results of fits performed on the data from RBC/UKQCD collaboration on pseudoscalar masses, 
decay constants and topological susceptibility \30f . The first fit lets F^^ vary freely. The other columns either a 
determination of the lattice spacings using the fl mass or 0{a^) discretisation effects for the decay constants. 
The star superscript indicates a quantity set to a fixed value (no uncertainty). These results should be compared 
with fits A2 and A5 in table 1. 19 



The results are given in Tab. [2] including finite-volume effects |j. First we consider fits B2, 3, 4 
without the topological susceptibility, including or not discretisation errors and lattice spacing de- 
termined from Mfi. Discretisation errors remain generally small (at most 5 %) apart for F^ in fit BA 
(20% effect, in good agreement with the results of ref. [30]), but lead to enlarged uncertainties on the 
other parameters. We notice that these discretisation effects are compatible with zero within error 
bars, which explains that the fits A2 and B3 (differing only through the effect of discretisation errors) 
yields very similar results. The (re) determination of lattice spacing through Mq performed in B2 
and B4: has a much more significant impact, as it tends to decrease the lattice spacings significantly 
(20-30%) as well as the value of the physical quark masses, and to increase Fj^/Fk noticeably. In 
these fits, the simulated quark masses stand much further away from the physical value than quoted 
in ref. [30j, implying that HO remainders at the simulated quark masses (scaling generically as p^ 
with respect to the HO at the physical point) play a significant role in the chiral expansion of observ- 
ables (up to 40% for the heaviest m^). The large error bars for dimensionful quantities is a refiection 
of the uncertainty on the lattice spacings determined from Mq. The significant difference between 
the value obtained the lattice spacings obtained here and in ref. [30] suggests a reanalysis the lattice 
data following the same procedure as the RBC/UKQCD collaboration (with the determination of 
the scaling trajectories and an iterative determination of masses and lattice spacings) but using the 
Re^PT formulae instead of analytic or standard xPT ones. Such a study would be very welcome to 
establish the effect observed in our fits, but is clearly beyond the scope of the present article. 

Finally, we include the topological susceptibility among the inputs in fit B5. As in Sec. [H the 
role of this input in the fit is marginal compared to the other data, and the outcome of fits B4 and 
B5 is very similar. Indeed, the Xmin obtained from fit B5 gets similar contributions from F,^, F^^M-^, 
FkMk and Mq, whereas the contribution from the topological susceptibility is small. The rest of the 
analysis is unchanged, with a competition between LO and NLO contributions for the observables of 
interest: 



^poie 

A. 



0.60 -F 0.25 -F 0.15, F'^Ml 

0.38 -h 0.52 -F 0.10, FJ^M]^ 

0.33 + 0.58 + 0.09 , F^M^ 
0.50 -F 0.35 -F 0.15. 



0.55 -h 0.30 -F 0.15, 
0.37 + 0.52 + 0.11, 
0.35 + 0.66-0.01, ^ ' 



The competition between LO and NLO terms of the chiral expansions observed in the previous fits 
remains though a bit less pronounced: 1^(3) is closer to one, and the value of Lg is closer to zero, 
but there is still an enhancement of NLO contributions to observables. Indeed, this enhancement 
is parametrised by msBoAL^ and itIsBqALq rather than L4 and Lg themselves, as can be seen 
for instance in eqs. ((IOl)-([II|) (with the typical values AL4 = L\{Mp) + 0.50 • 10~^ and ALg = 
Ll{Mp) + 0.25 • 10-3 for r = 25 and Y{2,) = 1) [UlIIS]. 

Our fit Bb features a fairly good Xmin/-^ compared to the results in ref. [30j, taking into account 
that our fit reproduces the physical values of the pseudoscalar masses and F^ by construction. We 
agree with ref. ^^ on the size of discretisation errors, but obtain different results concerning the 
determination of the lattice spacings from Mq. The main difference stems from our use of RexPT 
in the fit, as can be illustrated by performing the same fit as -B5 but constraining -^^(3) and Z{'i) to 
remain between 0.9 and 1.05. This last constraint mimics the usual assumption made in three-flavour 
xPT that both F^ and M,r are nearly saturated by their LO term. The minimum of the fit may 
look satisfying with r = 26.1, X(3) = 0.9, Y{2>) = 0.97, Z(3) = 0.93, Fk/F^ = 1.18 and values of 
the lattice spacings compatible at the 5% level with ref. [HC], but the value at the minimum is rather 
awful with Xmin ~ 151/11. This shows clearly that allowing for a numerical competition between LO 
and NLO is mandatory to reach decent fits to the whole set of data considered here. 

The results from fit Bh are also illustrated in the upper part of Fig. [5l The change in the value of 



^Contrary to Sec. 3] including or not finite-volume effects affects the value of Xmin only marginally, because the 
small changes due to finite- volume effects can be mimicked at the level of the fit by a slight shift in the values of the 
lattice spacings. This was also observed in the fits of the RBC/UKQCD collaboration on their own results |30) . 
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Figure 6: The pion decay constant as a function of the light-quark mass, according to the best fit in the last 
column of Tab.]^ (solid), with the same inputs but setting Y{3) = 1 and/or Z{3) — 1 (dashed and dashed-dotted 
curves), and with the same inputs but setting r = 35 (dotted curve). The physical point according to the fit 
result is indicated with a square. 

the lattice spacings rescales all lattice data points, which fall on the dashed curves corresponding to 
the best RexPT fit including finite- volume effects (responsible for most of the curvature at small m) 
more easily than in Sec. HI The values of p = rhs/ms and Y{3) are larger than in fit A5: from eq. (|10p 
(which has a similar structure to the Re^PT expression for F^^), we see that NLO contributions are 
proportional to the product pY{3) and tend thus to drive the (dashed) curves corresponding to the two 
volumes further apart than in the previous section. We also display the (solid) curve corresponding 
to the physical value of mg in an infinite volume, which by construction passes through the physical 
point Ftt = 92.2 MeV (this point corresponds to a slightly different value of m from Sec. |3|). The 
curvature is then essentially tied to the value at the origin, i.e., Z(^2): a lower value of Z(^2) will yield 
a steeper increase of Ft^ when rh increases. The dependence of Z(^2) on the various parameters can 
be read from the three-fiavour expansion of F^ in the chiral limit ruu = rrid = |31] : 



Z(2) 



r + 2 



[1 - riir)] + 



r + 2 



Z(3) 



r + 2 



y(3)5i + . . . 



(44) 



where gi denotes a small positive combination of chiral logarithms (around 7% near the physical 
point) and the ellipsis indicate HO remainders. Therefore, a dependence with a stronger curvature 
around the physical point and a flatter behavior above can be achieved by taking a larger value 1^(3) 
or a smaller value of Z{3), as illustrated in Fig. [6] and observed in fit B5 compared to A5. As an 
illustration, we indicate on the same figure the curves obtained with the same inputs as the best 
values for fit B5 of table [21 but setting Y{3) = 1 and/or Z(^3) = 1. The case of a complete saturation 
of the chiral series for F^ and F^M^ by their LO contribution (^(3) = Z{3) = 1) yields a higher 
Z{2) and a fiatter curve than our best fit. Finally, we notice that an increase of r (at fixed p) yields 
a slight increase of Z{2), but more importantly an upward shift of the physical value of rh, so that 
the corresponding line remains above the best-fit curve over a larger range of rh. The combination 
of these effects allows our formulae to reproduce a diversity of behaviours for the dependence of -F,^ 
on m, including the one exhibited by the RBC/UKQCD data. 

Since the results presented in table [1] do not include effects related to lattice spacings and since 
the determination of the lattice spacings involves assumptions on the form of chiral extrapolation. 
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we consider the most complete fit B5 in table [2] as the final result of our analysis, featuring a 
satisfying Xmin/-^- Interestingly, these results show a good compatibility with what was obtained 
in ref. [47j for the parameters of the LO chiral Lagrangian according to the results of the PACS-CS 
collaboration |40j . The value of the condensate in the Nf = 3 chiral limit is 

(S(3;2GeV))^/3 = 243±12MeV, (45) 

on the lower side of was obtained by the RBC/UKQCD collaboration (i.e., 256 it 6 MeV) [30], while 
the condensate in the two-flavour limit is 

(S(2;2GeV))^/3 = 285±8MeV, S(2; 2GeV)/S(3; 2GeV) = 1.51 ± 0.11 , (46) 

which illustrates the paramagnetic suppression of the Nj = 3 condensate with respect to the A'^^ = 2 
one (a similar statement holds for the pseudoscalar decay constant). 

6 Conclusion 

We have studied the topological susceptibility, a very interesting quantity related to chiral properties 
of QCD vacuum, and in particular to the quark condensate. This has led lattice collaborations 
to use this observable to determine the quark condensate, in addition to studies of the spectrum 
and dynamics of the light pseudoscalar mesons. However, it is important to assess higher-order 
corrections to the deceivingly simple connection between the topological susceptibility and the three- 
flavour quark condensate at leading order. Following our recent work on potential issues in three- 
flavour chiral extrapolations of lattice data [U], we have reassessed the information that can be 
extracted from this quantity allowing for a significant paramagnetic suppression of the Nf = 3 
quark condensate using the Resummed xPT (Re^PT) framework. We have noticed that for lattice 
simulations close to the physical situation (with a significant mass hierarchy between the dynamical 
strange and u, d quarks), the topological susceptibility would involve essentially the same combination 
of low-energy constants as the two-flavour quark condensate, and thus would not provide access to 
the three-flavour quark condensate as naively expected. In particular, using the leading-order three- 
flavour xPT formula would be particularly misleading if there is indeed a significant paramagnetic 
suppression of the three-flavour quark condensate. 

In order to escape this problem, two alternatives can be considered. A first possibility consists 
in performing further lattice simulations with hierarchy of quark masses different from the physical 
case, for instance with three dynamical quark masses of similar masses, or with a very significant 
strong isospin breaking. A second possibility relies on the combination of the topological suscepti- 
bility together with other sources of information on three-flavour chiral symmetry breaking, such as 
the spectrum of pseudoscalar mesons. Following this path, we have focused on RBC/UKQCD recent 
results [30] which provide data on pseudoscalar masses and decay constants as well as the topological 
susceptibility for two different volumes. Analysing the RBC/UKQCD unitary data points and in- 
cluding finite- volume effects, we have confirmed that these data do suggest such a suppression of the 
leading order Nf = 3 chiral order parameters (quark condensate and pseudoscalar decay constant) 
and the enhancement of next-to-leading order contributions related to the violation of the Zweig rule 
in the scalar sector (L4 and Lg contributions). The outcome of this first series of fits is mainly driven 
by the input from pseudoscalar masses and decay constants, and the data on topological susceptibil- 
ity shows a good compatibility with the pattern of three-flavour chiral symmetry breaking already 
found in ref. [17j . In addition, we confirm the difficulties (though at a lesser degree) encountered 
by the RBC/UKQCD collaboration to accommodate the dependence of Ft^ on the light-quark mass 
given by their data and the physical value of i*V at the same time. 

This problem has led us to reconsider the procedure used to determine the lattice spacing using 
our expressions to describe the dependence of the pseudoscalar observables on the quark masses. We 
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were not able to follow the same procedure as RBC/UKQCD for the determination of the lattice 
spacings (based on the determination of scaling trajectories and the dependence of the il. baryon 
on quark masses). However, we performed a joint fit of pion, kaon and Q observables to fix the 
lattice spacings. We also considered leading-order discretisation errors that may affect the kaon and 
pion decay constants. We finally performed a fit combining these two effects and adding data on 
the topological susceptibility. We obtained thus our final results given by the fit B5 in table El We 
noticed a significant enhancement (20%-30%) of the inverse of lattice spacings compared to the values 
quoted in ref. [30], as well as small discretisation errors (5% or less, and compatible with zero). Like 
in the previous analysis, the data on the topological susceptibility play a marginal role, but show a 
good compatibility with the rest of the data, yielding a satisfying Xmin/-^- Since they include the 
largest sets of data and use Re^PT consistently for the fit of the data and the determination of the 
lattice spacings, the results of fit B5 should thus be considered as the actual outcome of our analysis. 

The impact of our analysis on the determination of the lattice spacings is remarkable, and it calls 
for a confirmation by reanalysing the lattice data following the same procedure as the RBC/UKQCD 
collaboration (with the determination of the scaling trajectories and an iterative determination of 
masses and lattice spacings) but using the Re^PT formulae. However such a cross-check would go 
beyond the scope of the present article and our abilities. Let us add hat the inclusion of these effects 
does not modify the emerging picture of three-flavour chiral symmetry breaking already obtained 
in ref. [I7j, with a significant competition between leading and next-to-leading orders in the chiral 
series for masses, decay constants and topological susceptibility. 

We have used the RBC/UKQCD data as a case study for chiral extrapolations on the lattice, 
including both the pseudoscalar spectrum and the topological susceptibility, which can (and should) 
be extended to the results of other lattice collaborations. Our present analysis suggests also the 
inclusion of additional observables, e.g., other topological quantities such as the topological coefficient 
C4 as well as baryon observables, in our framework [26]. But even before considering this extension, 
the generation of further data points at lower quark masses would help clarifying the issue of the 
determination of lattice spacing as well as confirming the pattern of three-flavour chiral symmetry 
breaking emerging from our studies, with significant contributions from strange sea quarks leading 
to a non-trivial structure of Nj = 3 chiral expansions. 
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A NLO low-energy constants 

As recalled in ref. [U] and outlined in Sec. [21 the exact mass and decay constant identities obtained 
within the Re^PT framework for F^, F"^, ^n^w ^^^ ^K^K ^^"^ ^^ inverted to reexpress NLO LECs 
in terms of LO parameters of the chiral Lagrangian, physical quantities, and HO remainders: 

y(3)AL4 = — 1— ^[l-r?(r)-Z(3)-e], (47) 



(r + 2)MV 

F^ 

_L2L| 

8 Mr 



y(3)AL5 = l^[vir) + e'], (48) 



I F3 



^^^^^^^ = I^(7T^iit[i-^('')-^(3)-^]' (^9) 



if! 



Y\3)ALs = Y^jf2Hr)+d']. (50) 
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A. 
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rus 


P 


q 


24 


1.73 


0.005 
0.010 


0.04 


1.4980 


1.4707 


0.0211 
0.0342 


0.1098 


0.0962 


1.1420 


0.1924 
0.3111 



32 2.28 0.004 0.03 1.527 1.510 0.0163 0.1056 0.0962 1.0976 0.1539 



0.006 
0.008 



0.0232 
0.0302 



Table 3: Parameters of the unitary lattice sets taken from ref JWII. 



0.2196 

0.2859 



with e(r) defined in eq. ()14p . d, d' , e and e' combinations of remainders associated with the chiral 
expansions oi tt, K masses and decay constants respectively and 



r]{r) 



F 



K 



1 VF2 



1 



(51) 



ALj = L^ifJ-) — Li{jji) is independent of the renormalisation scale ^ and combine the (renormalized 
and quark-mass independent) constants -^4,5,6,8 together with chiral logarithms: 



o 
32^2^5(^) = \ 

327r^Le{fi) = 4 
16 



32^2^8(/^) 



o 2 

1. Mk 



1 



A^" 



log 



o 2 

Mk 
^2 



(r- l)(r + 2) 
o 2- 



o2 o2 

(4r + 1) log ^ + (2r + 1) log -^ 



+ 21og^ 



+ 



1 



o 2 



o2 

^ Mk 2, M„ 
log ^ + - log -^ 
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1 



o 2 

3 log — ^ + 5 log 
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Mk 
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16(r + 2)(r- 1) 
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o 2 



o2 

log ^ + - log -f 



+ 



16(r - 1) 



310^4 + 10^4 

M^ Mk. 

M^ Mk. 



(52) 
(53) 
(54) 
(55) 



where Mp are the LO contributions to the pseudoscalar masses, see eq. ([TH]) . 



B Lattice inputs 

We take our data points for pseudoscalar decay constants and masses as well as from topological 
susceptibility from the recent work of the RBC/UKQCD collaboration [30]. They considered 2+1 
dynamical flavours of domain wall fermions for two different lattice volumes 24^ x 64 x 16 and 
32^ X 64 X 16 (where the 16 corresponds to the extent of the fifth dimension inherent in the domain- 
wall fermion formulation of QCD). We consider only unitary sets where the masses of the sea and 
valence quarks are identical, with parameters recalled in Tab. [3] and observables in Tab. HI 

We denote Ag = a{rhq — rrires) the combination corresponding to bare masses (before addition of 
the residual mass rures and the conversion into the MS-bar scheme by a multiplication by Z™). We 
give a~^ and quark masses in units of GeV, x in units of 10"'' GeV , Fp in units of 10~^ GeV , FpMp 
in units of 10^^ GeV^. The dimensionful quantities have been converted from the lattice results by 
multiplying by the appropriate power of the lattice spacing, assuming for the latter the values quoted 
in the table. When the lattice spacings are allowed to vary and included in the parameters of the 
fits in Sec. [Sj these quantities are naturally rescaled by the appropriate power of the relevant lattice 
spacing. 
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V q 


X ■ 10-^ 


F^ 


F'k 


F^M'. 


F^Mi 


24 


1.1420 0.1924 
0.3111 


2.90(14) 

4.58(25) 


10.41(13) 
12.26(14) 


13.47(14) 
14.85(18) 


1.129(16) 
2.152(25) 


4.459(50) 
5.471(68) 


32 


1.0976 0.1539 
0.2196 
0.2859 


1.48(14)] 

3.22(55) 

3.63(72) 


9.799(96) 
10.65(10) 
11.53(10) 


12.99(11) 
13.51(11) 
14.16(12) 


0.820(10) 
1.265(13) 
1.788(18) 


3.988(36) 
4.390(38) 
4.895(42) 



Table 4: Lattice data for the pseudoscalar masses and decay constants as well as topological susceptibility 
taken from ref. 130f . 
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a-^ Aq 


As aMn 


24 


1.73 0.005 
0.010 


0.04 1.013(3) 
1.028(4) 


32 


2.28 0.004 
0.006 
0.008 


0.03 0.760(2) 
0.765(2) 
0.766(3) 



Table 5: Lattice data for the fl baryon taken from ref. 130]. 

In Tab. HI we do not include uncertainties coming from the determination of the lattice spacings, 
as these uncertainties would be completely correlated. Moreover, as explained at the end of Sec. HI 
we do not include the value of the topological susceptibility at a lighter quark mass given in ref. [30] , 
since it is likely to be affected by large systematics of unknown origin. In Sec. [5l we discuss the 
determination of the lattice spacings performed in ref. [30] using the mass of the il baryon, gathered 
in Tab. U 

A last comment is in order concerning the determination of the topological susceptibility in 
ref. [28], based on gauge configurations for a smaller volume (16)^ x32 in refs. |36p57j . In principle, this 
work could provide valuable additional information, but we have not been able to obtain consistent 
fits of the masses and decay constants of pions and kaons with the three ensembles (16)^ x 32, 
(24)^ X 64, (32)^ X 64, leading us to suspect an underestimation of the errors attached to the data for 
(16)'^ X 32. We have thus decided to keep only data for (24)^ x 64, (32)^ x 64, which were obtained 
for larger volumes and lighter quark masses, and thus less likely to be affected by sizable systematics. 

C Finite-volume effects 



As discussed in refs. |191l52f[55] and recalled in Sec. 13. 4^ for simulations where the time direction 
is much larger than the spatial ones, the finite-volume effects amount at NLO to a modification of 
the chiral (tadpole) logarithms. One can make contact with the infinite-volume integral (and the 
corresponding tadpole logarithm) in the following way, see eq. (j35p : 

o ^ o -^ 

ap Mp ^ Mp ° 2 

We have already analysed this expression in ref. |45j . ^1/2 contains a divergence of the form 

1/{L'^\J Mp) in the chiral limit. However, one expects ap to diverge like l/(L^Mp) due to the 
propagation of pion zero modes, and in ref. [45] , the alternative definition: 



Sp = ^-f|log^=^V2(^,Mp). (56) 



■=p 



Ci/2{L,M 



o 2 1 



P) 



L^^ Mp. 



^im;' <"' 
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was proposed to reduce the size of HO corrections to the chiral expansions (the term in square 
brackets is analytic in quark masses) and it was used to determine the contribution of finite-volume 
effects in Resummed xPT for pseudoscalar masses and decay constants. 

In the present article, we define the HO remainders of the chiral expansions in the limit L — )• c«, 
and we assumed that even at finite volume, the variation of these remainders with respect to the 
volume would be very small, considering that the simulations are performed in the p-regime where 
2mBQL^ ^ 1 |58H6U| . This is the reason for not including a volume-dependence on our expressions 
for the HO remainders. In a similar spirit, the choice between eqs. psp and ()57p amounts to a 
redefinition of what would be considered as HO terms, for instance: 



1 



4^22,3 



1 1 






+ 2 



1 1 



■ Mk. 



(58) 



inducing only a small numerical correction of the acceptable range of variation for these remainders 
in the p-regime. This would not be the case in the e-regime, where ImBoL'^ <^ 1, which deserves a 
separate study [26] and where the distinction between eqs. ([35]) and ([57]) could induce large differences 
in the acceptable range for HO remainders. 

Indeed, we have checked on the fits considered in the present paper that the outcome of finite- 
volume corrections according to eqs. (|35p and (|57p led to very similar results if we allow reasonably 
large ranges of variation for the HO remainders. For convenience, we quote only the results obtained 
with eq. ([35p . using for HO remainders the dimensional estimates described in ref. 
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